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Abstract 

A novel exactly solvable Schrodinger equation with a position-dependent mass 
(PDM) describing a Coulomb problem in D dimensions is obtained by extending 
the known duality relating the quantum d-dimensional oscillator and D-dimensional 
Coulomb problems in Euclidean spaces for D = (d + 2)/2. As an intermediate 
step, a mapping between a quantum d-dimensional nonlinear oscillator of Mathews- 
Lakshmanan type (or oscillator in a space of constant curvature) and a quantum 
D-dimensional Coulomb-like problem in a space of nonconstant curvature is derived. 
It is hnally reinterpreted in a PDM background. 
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1 Introduction 


In recent years, Schrodinger equations with a position-dependent mass (PDM) have at¬ 
tracted a lot of attention because PDM is of utmost relevance in a wide variety of physical 
situations, such as in energy density many-body problems, in electronic properties of semi¬ 
conductors and quantum dots, in quantum liquids, ^He clusters, and metal clusters (see, 
e.g., HI Eli for a list of references). As in constant-mass cases, exact solutions play an 
important role because they may provide both a conceptual understanding of physical phe¬ 
nomena and a testing ground for some approximation schemes. As a consequence, several 
methods have been tried to derive them. 

It is well known that Schrodinger equations in a PDM background may be alternatively 
interpreted as Schrodinger equations in curved spaces. On the other hand, a nonlinear 
oscillator in a space of constant curvature i El El [7], which has arisen as a quantization mu 
of the celebrated Mathews and Lakshmanan one-dimensional classical nonlinear oscillator 
[To] , has attracted much interest. The aim of this communication is to show that it may 
serve as a starting point for building new exactly solvable PDM Schrodinger equations. 

A clue to this construction is the relationship between the d-dimensional harmonic 
oscillator problem and the D-dimensional Coulomb one in Euclidean spaces for specihc 
values of the pair {d,D), resulting from coupling-constant metamorphosis |TT], Stackel 
transform |T5], regularization of the Coulomb problem [121 1131 [TB], supersymmetry [IB] , 
or duality tranformation im. Since the latter is a very simple and powerful method for 
generating new exactly solvable potentials (see, e.g., we choose to apply it in the 

present case. We plan to extend it to curved spaces, then to reinterprete the relation so 
obtained in a PDM background. 


2 Going from the oscillator to the Coulomb problem 
in Euclidean spaces 


As well known , the radial equation for the quantum d-dimensional oscillator in an Euclidean 
space reads 

f (f d — ld l{l + d—2) 2 2 ^ r, r, N 

- - - 5 - Oj^r^ + 2E ] R{r) = 0, 0 < r <+oo, (1) 

\dr^ r dr J 
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in units wherein h = m = 1. Here u and I denote the oscillator frequency and the an¬ 
gular momentum quantum number, respectively. It has an inhnite number of bound-state 
solutions, which, up to some normalization factor, can be written as 


oc r'e 


= 0 , 1 , 2 ,..., 


( 2 ) 


in terms of Laguerre polynomials, and correspond to bound-state energies 


d 


En = oj(n+-], n = 2nr + I- 


(3) 


The set of functions Rn^^i{r), corresponding to a given I value, is orthogonal on the half-line 
(0, -|-oo) with respect to the measure djj, = r^~^dr. 

On setting r = \fR and Rn^^i{r{R)) = Snr,i{R), we arrive at the differential equation 


_ I _ /(/ + d 2) ^ _ 1, ,2 ) c/'p'i _ Q 

^ 2RdR 4i?2 2R A ' ^ ’ 


(4) 


which can be interpreted as the radial equation for the quantum D-dimensional Coulomb 
problem, 

provided we set 


D=i(d+2), i = Q = \E, e = -\u]'‘. 


( 6 ) 


This shows that the roles of the coupling constant and the energy eigenvalue have been 
exchanged. On combining eqs. ([3]) and ([6]), we can rewrite Q as 

D - 1' 


Q = oj [u + 


where u = rir + L, 


(7) 


yielding oj = Q/ {u + Inserting this expression in that of the Coulomb bound-state 

energies given in ([6]), we get 






2(2z/ + D- 1)2’ 

with corresponding wavefunctions obtained from ([2]) in the form 


( 8 ) 




( 9 ) 
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The set of functions Snr,L{R), corresponding to a given L value, is orthogonal on the half¬ 
line ( 0 , -|-oo) with respect to the measure dfi = R^~^dR. 

Equations dH]) and (jHl) coincide with well-known results for the Coulomb problem in a D- 
dimensional space. It is worth stressing that the bound-state spectrum of the D-dimensional 
Coulomb problem is related to half the spectrum of the d-dimensional oscillator (namely 
that of even angular momentum states) for some even integer d, dehned in (IBl) . 


3 Going from a nonlinear oscillator to a Coulomb-like 
problem 

The one-dimensional classical nonlinear oscillator, first considered by Mathews and Laksh- 
manan [10], can be described in terms of a Hamiltonian 




( 10 ) 


where a plays the role of u in the standard oscillator and the nonlinearity parameter A 7 ^ 0 
enters both the potential energy term and the kinetic energy one. According to whether 
A > 0 or A < 0, the range of the coordinate x is (—cx), -|-cx)) or (—1/V^fAf, +1/ -^/jAf). 

The quantum version of H has been obtained [HIE] by replacing \/l + Ax^p by the 
operator —i\/l -|- Ax^ d/dx, yielding 


H = — 
2 


d 2 


d 




+ 




2(1 + Ax2)’ 


( 11 ) 


dx^ ' dx 

which is formally self-adjoint with respect to the measure djj, = {1 + Xx‘^)~^Rdx. Such a 
Hamiltonian is exactly solvable for a A-dependent potential parameter cR = /3(/5-fA). From 
now on, we will assume that o? is defined in this way. 

A d-dimensional generalization of the classical Hamiltonian ffTOj) has been proposed 
in such a way that the resulting Hamiltonian 


E = - 
2 




XiPt 


+ 


(1 + Ar2)^p2_^^j2 


i<j 


(3{I3 + A)r2 

2(1 + Ar2) 

^ 2(1 + Ar2) 


( 12 ) 
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keeps the maximal superintegrability property of the standard d-dimensional oscillator. 
In fir2|) . all snmmations rnn over i, j = 1, 2, ..., d, Jij = XiPj — XjPi denotes an angnlar 
momentnm component, and xf with r running on ( 0 , +oo) or ( 0 , 1 /according 

to whether A > 0 or A < 0. Furthermore, H may be interpreted as describing a harmonic 
oscillator in a space of constant curvature k, = —A. 

The quantization of flT^ in two HE] and three | 6 ] dimensions has been studied, but it 
can be easily extended to d dimensions. On replacing \/l + Ar^p^ and Jij by the operators 
—i\/l + Ar^ d/dxi and Jij = —i{xid/dxj — Xjd/dxi), respectively, we arrive at 


H = 


(1 + Ar^)A + Xr-^ + XJ^ 
or 


+ 


(3i(3 + xy 

2(1 + Ar 2 ) 


(1 + Ar^)^ + (d - 1 + dXr^)--^ - ^ 

CjT' ^ I CjT' 7^ ^ 


+ 


( 3{(3 + A)r2 

2(1 + Ar2)’ 


(13) 


with .P = Ylii<j Pj denoting the Laplacian in a d-dimensional Euclidean space. 

The corresponding Schrodinger equation is separable in hyperspherical coordinates and 
gives rise to the radial equation 


d^ 


(1 + Ar^)— + (d - 1 + dXr‘^)Jy 

dr^ rdr 


Id /(/ + d- 2 ) PP + X)r 


l + Ar^ +2Ofl(r)=0, (14) 


where P has been replaced by its eigenvalues /(/ + d—2), / = 0, 1, 2, .... The differential op¬ 
erator in (1T4|1 is formally self-adjoint with respect to the measure d/i = (1 + Ar^)“^/^r‘^“^dr. 
For d = 2, eq. ffTT|) reduces to eq. (29) of j?]. Its solutions can be easily obtained by 
extending the d = 2 approach to general d values and are given by 

Rn^r) oc r\l + + 2Xr^), n, = 0,1, 2,..., (15) 

in terms of Jacobi polynomials, with corresponding energy eigenvalues 


En = (3 [ n + ^ j — ^Ti{n + d — 1 ), n = 2nr + I- 


(16) 


The range of n values in flTbll is determined from the normalizability of the radial wave- 
functions on the interval (0, +cx)) for A > 0 or (0,1/ VTAf) for A < 0 with respect to the 
measure d/i = (1 + Xr'^)~^^‘^r'^~^dr. It is given by 

0 , 1 , 2 ,... ifA< 0 , 


n = 


0,l,2,...,7i 


— — +hl + 77, < — — ——1 if A > 0 

max 5 ^ 2 — '^hnax ^ ^ 2 j-j- ^ yj> 


(17) 
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It is worth observing that in the limit where /9/|A| goes to inhnity, the wavefunctions flT^ 
go over to ([2]) (with uj replaced by /3), dne to a limit relation between Jacobi and Laguerre 
polynomials [2U] . 

Let ns now perform the same transformation r = y/R and Rn^^i{r{R)) = Sn^^i{R) as in 
Sect. 2. This yields the differential eqnation 

+ ^(1 + + (d + l)Afi| A _ 

+ T|2B _ Xl{l + d-2)]- l/3(/3 + A) + iA£;)s(fl) = 0, (18) 

which can be rewritten as 


(1 + XRf 




+ 


D 


Q 


dR^ ' R 
S{R) = 0, 


(1 + XR) ( 1 + 


2D-1 \ d L{L + D-2) 

2D-2^^) Jr ^2 


(19) 


provided we set 


D = \{d + 2), L = \l, Q = \[E - 2XL{L + D - 2)1 8 =-\(3{(3 + X) + \XE. (20) 


It is straightforward to show that the differential operator in eq. (IT^ is formally self- 
adjoint on the interval (0, -|-cx)) for A > 0 or (0,1/|A|) for A < 0 with respect to the measnre 
djj, = {1 + XR)~^^^R^~^dR, corresponding to a space of nonconstant cnrvature. In snch 
a space, the potential —Q/R may not be interpreted as a Coulomb potential, since the 
latter, obtained as a solution of Laplace equation, assumes a more complicated form. We 
will therefore refer to it in this section as a Coulomb-like potential. 

As revealed by eq. fl20|) . the exchange of the roles of the coupling constant and the 
energy eigenvalue also looks less strict since L{L + D — 2) and E make their appearance 
in Q and 8, respectively. We can, however, proceed as in Sec. 2 and combine eqs. (IT6|) and 
fl20|) to write Q as 

Q = 13 H- 


- 1 


-A 


u \ u D — 2 L(L D — 2) 


u = rir + L, (21) 
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yielding {3 = {y + ^ \^Q + \ \v {v + D — + L{L + D — 2)] }. Inserting this expres¬ 

sion in that of S given in (1201) . we obtain after a straightforward calculation 


2{2iy + D - 1)2 

X / Q -|- A 


Q + A 


—z/ 1 L(^L D — 2) 


-{u + D - 1) [ u + D - -] + L{L + D - 2) 


( 22 ) 


It is worth observing here that, in contrast with the conservation of accidental degeneracies 
that occurs when going from the standard oscillator to the nonlinear one (see eqs. ([3]) 
and fITBD l and which is related to the maximal superintegrability property conservation, 
nothing similar happens in the Coulomb-like case since eq. (IE]) is replaced by fl22l) . As a 
matter of fact, this is an important aspect of our method. Using a more complicated A- 
dependent transformation might have maintained the maximal superintegrability property 
and, consequently, the accidental degeneracies, but it would not have given rise to any new 
result when going to a PDM picture. 

From the radial wavefunctions (IT5ll , we also get 


K R‘-(.l + Afi)"P<y>(l + 2\R), 


(23) 


where 


p = 2L + D - 2, 
1 


a = 


T = 


(^+ 

1 


Q + A 


u^ + iD- l)z/ + -{D -l) + L{L + D-2) 


(24) 


Q + A 


u ^ u D — — ) -1- L[L D — 2) 


X{2iy + D- 1) 

Bound-state wavefunctions, i.e., functions Snr,L{R) normalizable with respect to the mea¬ 
sure dp = {1 + XR)~^RR^~^dR, correspond to sets of quantum numbers {nr,L) satisfying 
the inequalites 


-|- (2L -\~ D — l)?7.j. -|- 2h^ -|- (2D — 3)17 -|- ^ 

nl + {2L + D — l)nr + L + -(^D — 1){2D — 3) < if A > 0. 


_g 

|A| 


if A < 0, 


(25) 


4' ^ A 

In both cases, there is only a hnite number of sets fulhlling these conditions and there is at 

least one {rir = 0, L = 0) provided Q > j{D — 1)|A| if A < 0 or Q > X(^D — 1){2D — 3)A if 
A > 0. 
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4 Reinterpretation in a PDM background 


As well known, one of the main difficulties of PDM problems comes from the non¬ 
commutativity of the momentum and mass operators, which can be coped with by using 
the von Roos approach [2T], wherein the kinetic energy operator is written as 
1 


E 


oxi oxi oxj axi 


(26) 


in terms of some ambiguity parameters rj, constrained by the condition .C + h + C = 

This form contains as special cases all the proposals that have been made in the literature 
and whose usefulness may depend on the physical problem in hand. To be more specific, we 
are going to consider here two special choices, namely the BenDaniel and Duke (BD) one 
(^ = C = 0 , r; = —1) [22] and the Mustafa and Mazharimousavi (MM) one (^ = C = ~l/4, 
V = -l/2) [3|. 

Considering first the d-dimensional nonlinear oscillator described by radial equation 
ffTTl) . it is straightforward to show that it can be reinterpreted as a d-dimensional nonlinear 
oscillator with a PDM m{r) = (1 -|- Ar^)“^, the BD and MM radial Schrodinger equations 
being 

did 


VAr) = 


Vi{r) - 2Eij Rnr.,iA) = 0) 

+ 


dr m{r) dr 

(' + ¥)(' + ¥) , /?(/?+A)r"-p 


1 -|- Ar^ 


2Ei = 2En - -d(d - 2)A, (27) 


and 


f-m ^/^(r)^m ^/“^(r) V 2 {r) - 2E^ Rnr,iA) = 0 ; 

\ dr dr J 


Voir) = 


r2 ^ _j_ 


2 E 2 = 2E,, 


(28) 


with Rn,,iA) = + XrA-^/^Rnr,iA) in both cases. 

Similarly, the D-dimensional Coulomb-like problem, characterized by radial equation 
(IT^ . becomes a D-dimensional (true) Coulomb problem with a PDM M{R) = (1 -|- XR)~‘^. 
The BD and MM radial Schrodinger equations now read 

did 


dR M(D) dR 


+ U{R)-2S,]SnMR)=0, 


(29) 










and 


+ U{R) - 2£^ Sn^R) = 0, (30) 

with the same potential 

rim (i + V)(i + ^) Q-K-D-1)(2C-5)A 

U(R) - -- - -, (31) 

and Snr,L{R) = + XR)~^^'^Snr,L{R), bnt different energy eigenvalues 

2£^ = 2£nr,L-U‘^D-l){2D-5)X^ 2 S 2 = 2Sn.,L - - S^X^ (32) 

As a final point, it is worth stressing that the PDM reinterpretation of Eqs. o and 
flT^ . which were associated with some complicated measures, has converted them in some 
equations in spaces with a simple measure dr or dR, which are directly applicable to those 
physical problems, mentioned in Sec. 1, wherein the PDM approach is relevant. 

5 Conclusion 

In this Letter, we have proved that the known duality relating the quantum d-dimensional 
oscillator problem to the quantum D-dimensional Coulomb one in Euclidean spaces, when¬ 
ever D = {d + 2)/2, can be easily extended to a quantum d-dimensional nonlinear oscillator 
of Mathews-Lakshmanan type (or harmonic oscillator in a constant curvature space). The 
result of the mapping is a quantum D-dimensional Coulomb-like problem in a space of 
nonconstant curvature. Going to an equivalent PDM description with respective masses 
m(r) = (1-|-Ar^)“^ and M{R) = {1 + XR)~‘^ leads to a duality between a nonlinear oscillator 
and a (true) Coulomb problem. If the existence of the former has been (at least implicitly) 
signalled in previous papers on low-dimensional problems [H El El IHl EH], that of the latter 
is a novel by-product of the simple transformation known in Euclidean spaces, which may 
lead to significant applications in those physical problems where the PDM concept appears. 

Another usefulness of the new exactly solvable model in multidimensional spaces, pro¬ 
posed in the present paper, may be in the context of the 1/N expansion approximation 
method [23l |24] . 
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Considering other duality transformations in curved spaces or in a PDM background 
than that connecting the oscillator-Coulomb pair would be a very interesting topic for future 
investigation. 
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